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ABSTRACT 

Despite decades of studying the Sun, the coronal heating problem remains unsolved. One funda- 
mental issue is that we do not know the spatial scale of the coronal heating mechanism. At a spatial 
resolution of 1000 km or more it is likely that most observations represent superpositions of multiple 
unresolved structures. In this letter, we use a combination of spectroscopic data from the Hinode 
EUV Imaging Spectrometer (EIS) and high resolution images from the Atmospheric Imaging Assem- 
bly (AIA) on the Solar Dynamics Observatory to determine the spatial scales of coronal loops. We 
use density measurements to construct multi-thread models of the observed loops and confirm these 
models using the higher spatial resolution imaging data. The results allow us to set constraints on 
the number of threads needed to reproduce a particular loop structure. We demonstrate that in sev- 
eral cases million degree loops are revealed to be single monolithic structures that are fully spatially 
resolved by current instruments. The majority of loops, however, must be composed of a number of 
finer, unresolved threads; but the models suggest that even for these loops the number of threads 
could be small, implying that they are also close to being resolved. These results challenge heating 
models of loops based on the reconnection of braided magnetic fields in the corona. 
Subject headings: Sun: corona — Sun: UV radiation — magnetic fields — methods: data analysis — 
techniques: spectroscopic 



1. INTRODUCTION 

While it is almost universally accepted that turbulent 
motions in the photosphere associated with magnetocon- 
vection provide the energy that heats the corona, there is 
no consens us as to whether this energy i s transferred by 
wave s (e.g. Ivan Ballegooijen et al.ll20lil : iMcIntosh et al.l 
l2011f h the dissipati on of magnet ic stresses through mag- 
netic reconnection (Parker 198 31), or injected direc tly by 
chromospheric jets (jDe Pontieu et all I2009L l2011f ). De- 
termining the details of the coronal heating mechanism 
would be greatly simplified if we could resolve individual 
coronal loops. The spatial resolution of instruments used 
to observe the corona, however, has generally been 1000 
km or more. The solar photosphere, in contrast, shows 
structure on scales of 100 km or less. We therefore cannot 
be sure that the properties we measure in the corona ac- 
tually correspond to those of individual loops. Arguably 
the most important question we need to answer is this: 
at what spatial scale do we resolve coronal loops? 

Previous work on this topic has shown that the 
temperature distribution in coronal loops is often 
narrow, suggesting that loops could be resolve d 
(jAschwanden fe Nightingale] 120051 : IWarren et all I2008D . 
These studies, however, either assumed that th ere is no 
substructure (Aschw anden fc Nightingale 2005) or sim- 
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ply computed a filling factor (Wa rren et al.ll2008h . The 
latter requires assumptions about the complex relation- 
ship between loop radius and measured width to crudely 
approximate the geometry. Here we use density measure- 
ments to infer the spatial scale of the emission and con- 
struct multi-thread models that can reproduce the cross- 
field intensity profiles observed with both EIS and higher 
spatial resolution AIA image data. With the geometry 
determined, we report the actual sizes of the strands. We 
generally find that a single thread implies loops that are 
much narrower than we observe, indicating that most 
loops are unresolved by current instrumentation. Our 
models show, however, that the observations can be re- 
produced with a relatively small number of threads sev- 
eral hundred km in width, and we do find a few cases 
of individual loops that are fully spatially resolved. This 
suggests that the ability to resolve coronal loops rou- 
tinely will be achieved by the next generation of coronal 
instruments. 

2. DATA REDUCTION AND ANALYSIS METHODS 

Understanding the observed intensity for optically thin 
coronal emission requires information about the physics 
of the atomic transitions being observed and the geome- 
try of the emitting structure. If we represent an observed 
loop as the superposition of a series of cylindrical flux 
tubes, the loop intensity is 



I = G(T,n)n 2 V/l 2 = G(T, n)n 2 Nitr 2 /I 



(1) 
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Fig. 1. — AIA 193 A filter image (upper panel) and EIS slit 
raster scans in Fe XI 188.216 A and Fe XII 195.119 A (center pan- 
els) showing a fully spatially resolved loop. This is loop #19 in 
Table [T] The traced loop segments are indicated. The best-fit EM 
distribution computed by the MCMC algorithm is shown in red 
in the lower panel. Random realizations of the solution are shown 
as dotted lines. The green curves are EM loci (observed intensity 
normalized by the G(T, n) function) for each line. Those lines that 
are not well correlated with Fe XII 195.119 A are indicated with 
dashed curves. The peak temperature, emission measure, electron 
density, and Gaussian width of the EM are given in the legend. 



where / is the observed intensity across a loop, G(T, n) 
is the normalized response of the plasma as a function 
of the electron temperature (T) and density (n), and V 



is the volume of the emitting structure within an instru- 
mental resolution element of area I 2 . The second expres- 
sion assumes that the emitting structures are a series of 
N identical loops of radius r. This equation allows us to 
derive information on the cross-field spatial scale of coro- 
nal loops, and theref ore to directly addres s our question. 

EIS on Hinode (jKosugi et all 120071 : ICulhane et"aH 
120071 ) is uniquely capable of determining the spectro- 
scopic diagnostic quantities required to perform this 
analysis. EIS routinely observes as many as 40 strong 
coronal emission lines in the wavelength ranges from 166- 
212 A and 246-292 A at approximately 22 m A spec- 
tral resolution. Several of these emission lines form 
density sensitive pairs which can be used to determine 
the electron density. They also cover a wide range 
of temperatures and their intensities allow us to com- 
pute the distribution of temperatures in a loop: the 
emission measure (EM) distribution. This methodol- 
ogy is well known and has already been applied to 
EIS observations of coronal loops by previous authors 
(iWarren et al.l 120081: iTVipathi et all 120091: iSchmelz et all 
l2Q10UDel Zanna et al.ll2011f ). 

AIA on the Solar Dynami cs Observatory (SDO, 

Pes nell et aHl20ll lL emen et al. 2011) observes the solar 
corona in 6 relatively broad wavelength ranges (1-11 A 
bandpasses). AIA has very high spatial resolution, 824 
km full width half maximum (FWHM) point spread func- 
tion ( PSF) compared w ith 1800 km for EIS (jGrigis et al.l 
l2012tlLang et al.ll2006[ ). but fewer plasma diagnostics. 

Here we collect a sample of 20 relatively isolated ac- 
tive region loops and use the AIA data to confirm the 
results we obtain with EIS. The EIS data have been pro- 
cessed for instrumental effects and converted to physi- 
cal units (erg/cm2/s/steradian) using the Solar Software 
(SSW) routine eis_prep. A specially designed observing 
sequence is used. This sequence uses the 1" slit to scan 
a field-of-view (FOV) of 240" by 512" in coarse 2" steps. 
The exposure time at each position is 60s. The duration 
of the scan is therefore about 2 hours and many diagnos- 
tic lines are included. Some of these are discussed below. 
The data are finally re-sampled to an arcsecond scale. 

AIA data were retrieved from the online cut-out ser- 
vice. These data are level 1.5 and have been corrected 
for flat fields and dark current. Bad pixels and cos- 
mic rays have been cleaned, and the data converted to 
DN/pixel/s. The data are internally coaligned by map- 
ping to the HMI (Helioseismic and Magnetic Imager) ro- 
tation and FOV center. Any residual offset is corrected 
by additional cross-correlation of images using the SSW 
routine align_cube_correl. The exposure time for each 
image is 1.7-4s, and we retrieved data covering the EIS 
FOV for plus / minus 2 hours around the time of the EIS 
loop observation. 

Affine transformation software developed for use with 
optical observations (which typically require higher accu- 
racy coalignment than possible for coronal images) was 
modified and used to obtain the relative scaling, rota- 
tion, displacement, and shear between the EIS Fe XII 
195.119 A raster images and the AIA 193 A filter image 
taken at the time the loop is observed. The EIS data are 
fitted to a double Gaussian function prior to this pro- 
cedure. The secondary Gaussian component is to take 
account of the weak blend at 195.18 A. Since the EIS 
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TABLE 1 



# 


Date/Time 




n 


T 


EM 


n 


TV 


r 


# 


Date/Time 


cr w 


n 


T 


EM 


n 


TV 


r 


1 


01/21 13:57 


1103 


9.68 


6.25 


27.6 


550 


6 


225 


11 


04/29 01:23 


1294 


9.12 


6.25 


27.3 


1450 


1 


1450 


2 


01/30 20:11 


859 


9.40 


6.30 


27.0 


690 


2 


404 


12 


05/06 13:55 


1208 


9.49 


6.20 


27.3 


650 


4 


330 


3 


03/18 09:34 


1468 


9.22 


6.20 


27.1 


935 


8 


335 


13 


06/14 00:52 


1222 


9.21 


6.25 


27.2 


1095 


5 


490 


4 


03/18 11:05 


1899 


9.13 


6.15 


27.5 


1980 


6 


815 


14 


07/02 03:32 


1388 


9.46 


6.20 


27.2 


620 


5 


280 


5 


04/15 00:41 


1526 


9.45 


6.20 


27.5 


850 


7 


325 


15 


07/02 03:53 


1186 


9.02 


6.20 


27.1 


1510 


1 


1510 


6 


04/15 01:41 


682 


9.46 


6.20 


26.7 


300 


1 


300 


16 


09/02 13:25 


1711 


9.41 


6.20 


27.7 


1040 


5 


470 


7 


04/21 12:02 


1426 


9.06 


6.30 


27.2 


1730 


3 


990 


17 


09/02 13:54 


1539 


9.28 


6.25 


26.9 


725 


6 


305 


8 


04/21 12:18 


1269 


9.20 


6.25 


27.5 


1255 


3 


730 


18 


09/17 18:57 


1859 


9.67 


6.25 


28.3 


1000 


5 


450 


9 


04/22 08:56 


1460 


9.21 


6.20 


27.3 


1220 


5 


550 


19 


10/14 22:55 


915 


9.06 


6.25 


26.6 


830 


1 


830 


10 


04/22 09:00 


1350 


9.25 


6.30 


27.5 


1210 


3 


700 


20 


10/14 23:51 


990 


9.58 


6.25 


27.3 


485 


3 


280 



The dates and times given indicate when EIS was rastering over the loop segment. All dates are in 2011. The parameter a w is the Gaussian 
width of the loop cross-field intensity profile in km measured in Fexn 195.119 A. The base-10 logarithm of the electron density (n) in cm , peak 
temperature (T) in Kelvin, and emission measure (EM) in cm -5 are given. r\ is the size of a single thread that matches the observed intensities 
in km. N is the minimum number of simulated threads and r is the maximum thread radius in km. 




Fig. 2. — Multi-thread simulations of the spatially resolved loop #19 shown in Figure [T] (upper panels) and the unresolved loop #14 
shown in Figure |4] (lower panels). The left column shows the simulated loop envelope (black circle) and the modeled threads (green). The 
observed cross-loop intensity profiles for EIS and AIA are shown in blue, and the simulated profiles are shown in red. A single thread 
completely fills the loop envelope and reproduces the intensity profiles for loop #19. Five 280 km threads are needed to reproduce the 
profiles for loop #14. 



pointing varies with the satellite orbit, the sampling in 
the slit rasters is not uniform across the FOV. Therefore, 
any inaccuracy in the AIA-to-EIS image registration is 
corrected visually as a final step. The EIS and AIA im- 
ages shown in the Figures in this Letter were treated 
with a Gaussian sharpening filter with a radius of 3 pix- 
els for presentation, but all analysis was performed on 



the unfiltered data. 

Suitable loops were first identified in the EIS raster 
data, and relatively isolated portions selected. Once 
identified, cross-loop intensity profiles perpendicular to 
the loop axis were automatically extracted and then 
averaged along the loop within t he selected segment 
following established methodology ([Warren et al.l 12008: 
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lAschwanden et al1l2008[ ). We then subtract a co-spatial 
background component by visually identifying the loop 
in the intensity profile and fitting a 1st order polyno- 
mial between two background pixels. This procedure is 
applied to the strong Fe XII 195.119 A line data, and 
the same background locations are used for all the other 
spectral lines. For this analysis, we use spectral lines 
from consecutive ionization stages of Fe VIII-Fe XVII 
that we have found to be reliable in previous EM analysis 
(jWarren fe Brooks! I2QQ9I : iBrooks et "all l2QQ9h . A Gaus- 
sian function is fitted to the background subtracted in- 
tensity profile to obtain the loop intensity and FWHM. 
The electron density is then calculated using the rec- 
ommended Fe XIII 202.044/203.826 ratio, and the EM 
analysis is made using all the other lines. The CHIANTI 
database is used for the atomic data includ ing the new 
ionization fracti ons (Q Dere et al. 1997, 2009). Following 
previous work ([Warren et al.l 120081 ) we only use those 
lines, the intensity profiles of which are highly correlated 
with that of Fe XII 195.119 A. The instrumental cali- 
bration uncertainty of 22% is used as the measurement 
error of these lines (|Lang et al.l [2006). The other lines 
have their intensities set to zero, and their uncertainties 
set to 20% of the background emission. Finally, for each 
loop we measured the FWHM of the cross-field intensity 
profile in both EIS and AIA. 

The EM calculation is performed using a Monte 
Carlo Markov Chain (MCMC) algorithm available in the 
PINT ofALE software package ([Kashvap fe Drake! 119981 
2000). 100 realizations of the solution are calculated. 
In addition, we separately compute another 100 Monte 
Carlo simulations to verify that the best-fit emission 
measure, thermal width, and peak temperature do not 
change significantly if we perturb the intensities ran- 
domly within the calibration uncertainty. As a consis- 
tency check, we also fit the data with isothermal and 
Gaussian emission measure functions. Finally, we en- 
sure that th e selected loops sat i sfy t he width criteria 
suggested by Aschwande n et all (j2QQ8h . but we do not 
adopt their background contrast criteria because we are 
interested in detecting cases where an EIS loop is com- 
posed of several AIA loops. Such cases would be missed 
if the procedure is optimized to detect single AIA loops. 

The multi-thread models are calculated at a resolution 
of 1 km using Equation [T] assuming an isothermal plasma 
for each thread and that the threads are of equal radii. 
The threads are evenly distributed within the loop en- 
velope. The thread intensity profiles are combined to 
produce a composite profile for the loop, which is con- 
volved with the instrument PSF and resampled to the 
instrument pixel scale. This profile can be compared 
to the observations. A broad parameter space of loop 
envelope and thread radia is explored to find the best 
match with observed EIS loop intensities and FWHM. 
Finally, a minimum number of threads compatible with 
both the EIS and AIA data is determined and their posi- 
tions slightly perturbed to obtain better agreement with 
the asymmetric intensity profiles. 

3. RESULTS AND DISCUSSION 

Representative images and calculated EM distributions 
for the loops in our sample are shown in Figs. [T] and 
IU The electron density, peak temperature, peak EM, 
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Fig. 3. — Ratios of simulated and observed intensities and 
FWHMs for both the EIS and AIA measurements for the entire 
loop sample. The simulations are able to reproduce the observed 
values to within 20% for most of the loops. 

and the Gaussian width of the loop cross-field intensity 
profile are presented in Table [TJ In all cases the com- 
puted EM distribution is a strongly peaked function of 
temperature. The narrow temperature distributions are 
consistent with p r evious results using spectroscopic data 
(jDel Zannall2QQ3t IWarren et al.l 120081) and mor e limited 
imaging data (jAschwanden fe Nightingale 2005), though 
the coverage of the imager observations was not sufficient 
to produce convincing temperature measurements. None 
of these studies were able to say anything definitive about 
possible substructure, but the very narrow temperature 
distributions measured for these loops suggest that they 
may consist of single monolithic structures. From Equa- 
tion [T] we see that we can use the EIS measurements of 
the density and temperature to infer the physical radius 
of the loop for the case of a single thread (7V=1), which 
is 7*1 in Table [TJ We emphasize that measuring the den- 
sity is critical for probing the internal structure of a loop 
since it allows us to isolate the geometrical factors that 
give rise to the observed emission. 

The application of this procedure shows that for a mi- 
nority of the cases that we have studied the EIS and AIA 
loop widths predicted from a single thread model are able 
to reproduce the observations. An example is loop #19 
which is shown in Figs. [TH21 This loop is compatible 
with a simulation of a single thread with a radius of 830 
km, and our emission measure analysis indicates that 
it is very nearly isothermal. Our results are the first to 
demonstrate unambiguously that some coronal loops can 
be resolved with current solar instrumentation . A pre- 
vious study (Aschwa nden fe N ightingale 2005) showed 
that the observed loop FWHMs are often larger than 
the instrumental PSF, as we also find here. This analy- 
sis, however, lacked density measurements and a model 
to deconvolve the physical size of the loops from the ob- 
served widths and instrumental PSF. As the authors of 
that paper pointed out, their analysis assumed a filling 
factor of unity, which is inconsistent with t he measured 
values for coronal loop s (Warren et al. 2008). In a follow- 
up paper, lAschwanden et all (2007) did model the sub- 
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resolution and considered broadening due to the PSF, 
but, without a density diagnostic, could not measure the 
volume of emitting plasma in coronal loops. 

For the majority of the loops in our sample, however, 
a single thread model produces simulated loops that are 
much narrower than what is observed. Such loops re- 
quire substructure to reproduce the measured FWHMs. 
For these loops we have determined the minimum num- 
ber of threads compatible with both the EIS and AIA 
data. The number of threads and the thread radius is 
given in Table HJ For most cases, the data are modeled 
to within the uncertainty of the instruments' calibration 
(20-25%). The ratio of the simulated and observed in- 
tensity and loop FWHMs for each loop is summarized 
in Fig. 03 Since the multi-thread model has been con- 
strained by the EIS observations the good agreement be- 
tween the EIS measurements and the model is expected. 
The model parameters, however, are not derived from the 
AIA data and this comparison serves as an independent 
confirmation. As a further consistency check, we have 
verified that there is good agreement (within 10%) be- 
tween the measured FWHM for Fe XII 195.119 A and 
the FWHM of the other lines. The models also still 
reproduce the intensities and FWHMs if the density is 
perturbed by the uncertainty in agreement between ob- 
served and calculated intensities (~ 35%). Significantly, 
although the strand radius may change by this amount, 
the number of threads remains the same and the size of 
the loop envelope does not change appreciably. 

An unresolved loop (#14) is shown in Figs. [2] and 
[H Here 5 threads with a maximum radius of 280 km 
can reproduce its observed properties and the emission 
measure analysis indicates that it has a finite but narrow 
thermal width, consistent with it being composed of a few 
threads. This loop is more typical of the majority of the 
sample. Our observations and models suggest, however, 
that the true size scale of even these coronal loops is not 
far from being resolved. A maximum of only 8 threads is 
needed to reproduce all of the loops in our sample, and 
the mean maximum thread radius of 480 km is not much 
smaller than the 720 km spatial pixels of EIS. Even the 
loops in our sample with the thinnest threads could be 
resolved with an instrumental spatial resolution of 200 
km. This seems an achievable goal and we suggest that 
these results could be a reference to justify the spatial 
resolution required for a next generation spectrometer. 
Such an instrument would be able to routinely resolve 
coronal loops for the first time. 

We stress that there is some ambiguity in the inter- 
pretation of these spatial scales since we do not know 
the maximum number of threads or the minimum thread 
radius. Loop #14 in the figures is an intriguing exam- 
ple, however, since the EIS intensity profile is reproduced 
with relatively few threads, that are then seen to sepa- 
rate out when observed by AIA. It is difficult to justify 
large numbers of threads if the profile is not smooth. 
Furthermore, if they have to be constrained to specific 
positions the profile would be dependent on the viewing 
angle. Observations from the STEREO A and B space- 
craft do not show any e vidence that loop prope rties are 
line-of-sight dependent ([Aschwanden et aLll2008[ ). 

Furthermore, our estimates of the loop width (r, n) 
are also supported by very high spatial resolution (100 
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Fig. 4. AIA 193 A filter image (upper panel) and EIS slit 
raster scans in Fe XI 188.216 Aand Fe XII 195.119 A (lower panels) 
showing an unresolved loop. This is loop #14 in Table ^ The 
traced loop segments are indicated. The EM analysis of this loop 
is shown in the lower panel. The description of this panel is the 
same as that of Figure [T] 



km) optical observations of "coronal rain", which are 
condensations that fo rm during the catastrophic cool- 
ing o f coronal loops ([Antolin fc Rouppe van der Voortl 
l2012f ). Such events may be rare, but they offer unique 
insights into the structure of coronal loops. These ob- 
servations yield a typical FWHM of 310 km, equivalent 
to a completely filled loop with a 215km radius when 



6 



Brooks et al. 



convolved with the instrumental PSF. They also show 
that neighboring cooling features occur nearly simulta- 
neously, implying a coherence to the evolution of separate 
threads that is consistent with the narrow EM distribu- 
tions we find. Such narrow distributions also suggest 
that the number of threads is small. An important av- 
enue of future research is to investigate the evolution of 
the resolved loops to see if they are consistent with sim- 
ple hydrodynamic mode ls, or are heated q uasi- steadily 
([Brooks fc Warrenl 120091 : iMulu-Moore et al.ll2011[ ). 

Our observations present a challenge to current views 
on the fine magnetic coupling of loops to the pho- 
tosphere, and heating models based on Parker's idea 
of nanoflare reconnection of braided magnetic fields in 
the corona (jParkerl Il983). Typical magnetic flux and 
field strengt hs suggest ~ 100 km radius photospheric 
flux tubes ([Priest et al.l 12002). and together with the 
very small spatial scales expected for magnetic recon- 
nection, suggest that coronal loops observed at 1000 km 
scales should be composed of many hundreds of inde- 
pendent threads in various stages of heating and cooling 
(Kli mchukl l2006f ) . This implies a broad distribution of 
temperatures at the spatial resolution achieved by EIS 



and AIA (Schmelz et al. 2011). Localized area expansion 
of the magnetic field could bridge t he gap between pho - 
tospheric and coronal spatial scales ([Warren et al. 2010), 
but there is a clear need to investigate the mapping of 
specific loops to the measured magnetic field. Finally, 
our observations of resolved single loops and EIS bun- 
dles reproduced by a few threads that split out into spa- 
tially separated individual structures when viewed by 
AIA, suggest that either magnetic reconnection occurs 
on much larger spatial scales than previously imagined, 
is energetically unimportant for loop heating, or the mag- 
netic field in the corona is generally topologically simple 
and not braided. 

We would like to thank the referee for their posi- 
tive comments and suggestions. This work was per- 
formed under contract with the Naval Research Labo- 
ratory and was funded by the NASA Hinode program. 
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ISAS/JAXA, with NAOJ as domestic partner and NASA 
and STFC (UK) as international partners, it is operated 
by these agencies in co-operation with ESA and NSC 
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